Estimation of the ratio of the pn^pn7T°7T° /pn^dn^TT^ cross sections 
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Abstract 

Evidence has recently been presented for the existence of a dibaryon of mass 2380 MeV/c 2 and width 70 MeV/c 2 , which decays 
C| strongly into the d7r°7r° channel [M. Bashkanov et al, Phys. Rev. Lett. 102 (2009) 052301; P. Adlarson et al., arXiv:1104.0123]. 

The decay rate of such a hypothesised dibaryon into the {pn}j=( ) ir°ir channel is estimated in a weakly model-dependent way by 
,-— J using final state interaction theory. It is shown that, if the resonance exists, it should then show up as strongly in this channel as 

in div n . The sum of the two decay modes would saturate most of the inelasticity predicted in the relevant partial waves in the 
£h 2380 MeV/c 2 region. 
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The experimental search for dibaryon resonances has 
been a long and generally painful story, with most claims 
being eventually disproved. However, the WASA collab- 
oration working at CELSIUS has recently produced evi- 
dence for a very significant peak in the total cross section 
for the quasi- free pn — > c?7r°7r° reaction [1]. The experi- 
ment was carried out using a deuterium target, with the 
final deuteron and the two pions from the pd — > pdir^ir 
reaction being detected and the proton in the final state 
reconstructed from energy-momentum conservation. Slow 
protons were then selected so that they could be treated as 
spectators that were assumed only to influence the reaction 
through the kinematics. Furthermore, because the events 
were fully constrained, the centre-of-mass (CM) energy W 
of the pn system could be determined on an event-by-event 
basis without any uncontrolled smearing arising from the 
deuteron Fermi motion. In this way the peak region could 
be scanned with a proton beam of fixed energy. 

After the move of the WASA facility to COSY-Julich, the 
experiment was repeated with the same basic apparatus but 
with significantly higher statistical precision, which allowed 
much finer divisions in the CM energy to be presented [2] . 
In order to ensure that only low spectator momenta were 
used in the subsequent analysis, data were in practice taken 
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at several proton beam energies and, although there is some 
uncertainty in the relative normalisation of these runs, the 
good overlap lends confidence that the spectator distribu- 
tion was correctly handled when deriving cross sections. A 
parameterisation [3] of the COSY results is shown in Fig. 1 
in terms of the excess energy Q = W — (nrid + 2m 7T o)c 2 . 

The pn — > dir°Tr° total cross section is dominated by a 
very strong peak at W « 2.38 GeV, i.e., Q w 234 MeV, and 
width r w 70 MeV, which the authors have speculated to be 
a signal for a dibaryon resonance. Since a 7r°7r° system can 
only have isospin / = or / = 2, such a resonance, as well 
as the dipion itself, must be isoscalar. The likely dibaryon 
states would therefore be J p = 1 + or 3 + , which couple to 
the initial ( 3 Si, 3 Di) or ( 3 Ds, 3 Gs) proton-neutron waves, 
respectively. 

The other important point to note is that the shape of 
the two-pion spectrum from the pn— »d7r°7r° reaction is far 
from that of phase space, with an enhancement at very low 
7r°7r° masses and also some excess of events at maximum 
ro^o^o . These are features that are typical of many two-pion 
production reactions in the isoscalar channel and go by the 
name of the ABC effect from their original observation in 
the pd — > 3 He(7T7r)° reaction [5]. 

The ABC effects are generally thought to be dynamical 
in origin, often associated with the production of two p- 
wave pions, such as through the excitation and decay of 
two A(1232) isobars [6]. However, it is hard to construct a 
model of this type that would give a very strong peak in the 
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Fig. 1. Total cross section for the pn^dn°ir° reaction as a function of 
the excess energy Q. The solid (blue) curve represents a parameter- 
isation [3] of the COSY- WAS A data [2]. The cross section predicted 
for pn^pnn°TT° on the basis of an L = 2 decay using the ratio of 
Fig. 2 is shown by the dot-dashed (red) curve and the sum of the two 
components by the higher dot-dashed (black) curve. The inelastic 
cross sections in the 3 D^ (magenta) and the 3 Ds + 3 G3 (turquoise) 
waves taken from the SAID current solution [4] and scaled by an 
isospin factor of 1/6 are shown as dashed curves. 

overall centre-of-mass energy, with a width that is only half 
that of a single A(1232). The most complete evaluation in 
such a model in fact gives T 200 MeV [7] . Even though 
the existence of a dibaryon might seem very improbable, 
the possibility has nevertheless to be taken seriously. With 
that in mind, one then has to ask if the dibaryon could and 
should be seen in other decay modes. It is the purpose of 
the present letter to evaluate the rate for the decay into 
the prnr°ir° channel in a weakly model-dependent way by 
using final state interaction theory. 

We have previously shown that at low proton-neutron 
energy E pn there is an intimate connection between the 
deuteron bound-state wave function ip a (r) and that of the 
pn 3 Si scattering state ip(q, r) [8]. Although it is well known 
that the continuation of the latter into the bound-state 
region does yield the deuteron pole, it is less obvious that 
the residue at the pole is uniquely fixed by the deuteron 
binding energy B. The connection is given by 

Mg,r)l 2 * . I* 2 , l^(r)| 2 , (1) 
a(a z + q z ) 

mB, with to the average nucleon mass and 
pn . The relation in Eq. (1) is exact at the pole 



where a 2 
q 2 = mE, 



q = —a but at short pn separations it remains a good 
approximation up to energies of at least 50 MeV [9] . 

The production of a meson in nucleon-nucleon collisions 
necessarily involves a large momentum transfer which is 
sensitive to short-range effects in the NN system. On the 
basis of Eq. (1) it is easy to see that for a short-range pro- 



duction operator the absolute squares of the matrix ele- 
ments leading to the dX and {pn} q X final states are re- 
lated by 



\M(pn -> {pn} q X) 
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a(q 2 + a 2 ) 



\M(pn^dX)\ 2 



(2) 

where {pn} q denotes a spin-triplet S*-wave proton-neutron 
state at relative momentum q. Deviations arising from the 
pn tensor force and the deuteron D-state are also likely to 
be small at short pn separations. 

The form of the relation in Eq. (2) depends upon the 
normalisation used. Here we have taken uniformly a boson 
normalisation, where the relativistic n-body phase space is 
given by 



d<£ n (p, P1 ,--- Pn )=5 4 p-$> n 



d 3 Pl 
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where (E^p^ is the four- momentum of one of the final 
particles with mass to^ and P is the total four- momentum. 
The differential cross sections then correspond to 
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where the flux factor F is the same for the different n-body 
final states. 

Using this approach, it is possible to describe quite 
accurately the ratios of the cross sections for single- 
meson production in pn — > {pn} q n/dn [10] and pp — > 
{pn} q ir + 1 dir + [9] close to threshold. In the following we 
carry out a similar analysis for two-pion production in the 
WASA experiments [1,2]. 

Whereas a 1 + dibaryon decay to dTT°TT° could involve 
purely 5-waves, that of a 3 + state necessarily requires L = 
2 waves. Because the observed two-pion spectrum is peaked 
in the s-wave region of low 7r°7r° masses, the angular mo- 
mentum must then correspond primarily to one between 
the deuteron and the dipion, where the relative momentum 
is k. This introduces a kinematic factor of k L in the near- 
threshold amplitude and this must be included in the mod- 
elling of the cross sections. We treat separately the cases 
of L = and L = 2 in order to estimate the ratios of the 
total cross sections: 



R L (Q) = a(pn -> {pn} I=0 Tr"ir {> )/a(pn -> dir ^ 



(5) 



If the non- relativistic approximation is made to the phase 
space of Eq. (3) and one neglects any dependence of the 
matrix elements upon the kinematic variables apart from 
the k L factor, it is possible to obtain estimates of the cross 
section ratio of Eq. (5) in closed form for both values of L: 
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where x = Q/B and B is the binding energy of the 
deuteron. 

The predictions of Eqs. (6) and (7) are shown in Fig. 2. 
The curve corresponding to L = lies significantly higher 
than that of L = 2, though some of this difference might be 
reduced a little if barrier penetration effects were included 
in the L = 2 case [11]. For the four-body final state some 
of the energy is taken by the pn internal degrees of freedom 
and this reduces the influence of the threshold k 2 factor in 
the L = 2 case. The curves increase with Q because of the 
larger phase-space volume there for a four-body final state. 
It is clear that at the indicated WASA peak position the 
pn^pnir°TT and pn— >d7r°7r° total cross sections should be 
of the same order of magnitude. 
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Fig. 2. Predicted ratio of the pn^rpniT n° /pn^dn°n° total cross 
sections as a function of the excess energy Q for both the L = 
and L = 2 possibilities. The dashed (red) curves result from the 
non-relativistic analytic formulae of Eqs. (6) and (7), where the tt°tt° 
spectrum is described by phase space. These predictions are reduced 
slightly to the dot-dashed (blue) curves when relativistic kinematics 
are used for the pions. The introduction of an ABC-like structure in 
the 7r°7r° effective-mass distribution increases the ratio to the solid 
(magenta) curves. The position of the centre of the WASA peak [1,2] 
is indicated by the (black) vertical dashed line. 

It is, however, clear that for Q in the 200-300 MeV range 
it is dangerous to treat the pions non-relativistically, even 
though many effects cancel in the cross section ratio. How- 
ever, the integration over a relativistic phase space does not 
lead to analytic forms analogous to those of Eqs. (6) and 
(7) . By retaining non-relativistic kinematics for the nucle- 
ons, the ratios can be expressed in terms of two-fold inte- 
grals which have to be evaluated numerically. This proce- 



dure reduces slightly the predictions of the non-relativistic 
approach, as shown in Fig. 2. The difference can be under- 
stood because, on average, the pions in the pnir°7r a final 
state are slightly less relativistic than those in d7r°7r°. 

In the relativistic evaluation thus far presented, the two- 
pion spectrum was still assumed to have a phase-space form 
whereas the experimental data show a strong tendency for 
a peaking at low tt°tt° invariant masses [2]. This is par- 
ticularly evident at the cross section maximum, but gets 
softer at higher beam energies. To investigate the effect, the 
7r°7r° spectrum at the "resonance" position was fitted and 
it was then assumed that this shape scaled with the avail- 
able excess energy Q. This assumption brings the problem 
closer to one of a three-body/two-body comparison [9,10] 
and increases the ratios to be equal or even above those of 
the non-relativistic formulae shown in Fig. 2. However, it 
is comforting to note that the evaluations of the ratios in 
the various approximations do not give significantly differ- 
ent results. The only major dependence appears to be on 
the value assumed for L in the decay. 

The predicted cross section ratio is multiplied by a pa- 
rameterisation of the WASA pn — » d7r°7r° total cross sec- 
tion [2] in order to get estimations for that of pn — > pnir°ir°. 
For this purpose no attempt was made to model the cross 
section peak in terms of a resonance and background be- 
cause events from the background are likely to yield a 
broadly similar four-body / three-body ratio. The results are 
shown separately for L = 2 and L = in Figs. 1 and 3, 
respectively. It is clear from the figures that a peak in the 
pn — > d7r°7r° cross section should also lead to one in pn — > 
pnir°ir at a marginally higher value of Q and very slightly 
broadened. 
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Fig. 3. Analogous curves to those in Fig. 1 but for an L = decay 
of the dibaryon. These are compared to the inelastic cross sections 
in the two J = 1 partial waves. 

Also shown in the figures are the sum of the three- and 



•3 



four-body cross sections. It is then of interest to see how 
these compare with the inelasticities predicted in a current 
partial wave analysis. It is important to note here that the 
pn — > dTT°Tr° represents only one sixth of a pure isoscalar 
1 = cross section and so the SAID predictions [4] for the 
allowed partial waves have been reduced by a factor of six 
before being plotted in Figs. 1 and 3. 

Although it is self-evident that the SAID predictions 
can only be as reliable as the experimental input, which is 
rather sparse in this energy range, the figures suggest that 
L = decay is not favoured because the estimated cross 
section significantly exceeds the sum of the SAID inelastic 
cross sections in the 3 Di and 3 Si partial waves [4]. This 
point is reinforced by the fact that some of this inelasticity 
will correspond to single-pion production in, for example, 
pn — > ppir~ . The same argument regarding the inconsis- 
tency with the SAID predictions cannot be made against 
the J p = 3 + possibility that is illustrated in Fig. 1. This 
preference is consistent with the 3 + deduced by the exper- 
imental group on the basis of information extracted from 
the angular distributions. 

Although the cross section ratio that we have evaluated 
has been presented in terms of a dibaryon resonance, the 
work here does not really rely on this hypothesis; it de- 
pends principally upon the short-range nature of the two- 
pion production operator. The curves in Figs. 1 and 3 have 
therefore a much wider range of applicability and the state- 
ment regarding the strong preference for the J p = 3 + waves 
dominating the cross section remains. On the other hand, 
if one did believe in the existence of such a dibaryon then 
our work could be interpreted as finding a relation between 
the dTT°TT° and pnir°ir° couplings to this state. 

We have only estimated one particular inelastic channel, 
viz. the pn7T 7r°, where the proton- neutron pair emerge in 
the spin-triplet S-wave. There are other potential channels 
for the decay of a J p = 3 + dibaryon, such as pptt~tt°, but 
their estimation falls outside the remit of our approach and 
would require a full dynamical model. Although one might 
question the use of the extrapolation theorem of Eq. (1) 
for Q in the 200-300 McV range, in fact the bulk of the 
energy is taken by the two light particles and the fraction 
where the pn excitation energy is large enough to jeopardise 
the approach is comparatively small. There are probably 
greater uncertainties associated with the assumptions of 
the dependence of the matrix element upon the possible 
kinematic variables. However the strong distortion of the 
7r°7r° invariant mass spectrum only changes predictions by 
typically 10%. The angular momentum in the decay seems 
to make the much bigger difference illustrated in Fig. 2. 



Independent of the origin of the peak in the pn — > dir°n 
total cross section identified by the WASA collabora- 
tion [1,2], the approach proposed here shows that there 
should be an analogous peaking in pn — > pnir°Tr° cross 
section of a similar strength. This channel could be in- 
vestigated using the existing WASA data, though the 
kinematics of this state are harder to identify due to the 
presence of the final neutron [12]. 
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